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We report results from a systematic strong-coupling expansion of a spin-i Heisenberg chain cou- 
pled to Einstein phonons. In the non-adiabatic regime (hil ~ J) this model is used to describe 
zero temperature properties of CuGeOs. The linked cluster expansion allows the determination of 
observables in the thermodynamic limit preserving the full lattice dynamics without a truncation 
of the phononic Hilbert space. In particular, the spin gap and the dispersion of the elementary 
triplet excitation are calculated up to 10 th order in a dimer expansion. The magnetic structure 
factor of the ground state is evaluated up to 6 th order. We show that the spin-phonon coupling 
leads to a renormalization of the elementary triplet dispersion. Surprisingly in the non-adiabatic 
regime a substantial renormalization of the spin gap only sets in at much larger couplings than those 
proposed for CuGeOs. The ground state magnetic correlations are found to be hardly effected by 
the spin-phonon coupling, but dominated by the frustrating magnetic interaction in the parameter 
regime relevant for CuGeOs. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.50.Ee, 71.38. +i, 63.20.Kr 



A renewed interest for the magnetic properties of one 
dimensional spin chains was created by the observation of 
a spin-Peierls (SP) transition in the inorganic compound 
CuGe03ti The SP phase of CuGeC>3 is characterized by 
an energy gap in the spin triplet excitation and a dimer- 
ization of the lattice along the chain direction. Experi- 
mentally, the nature of the SP transition was confirmed 
by inelastic neutron scatterinSr-(iNS|)-,rSusceptibility and 
X-ray diffraction experimentso'tfffl'ErLlQ. 

The SP transition is usually driven by the coupling of 
spins to the lattice. Previous theoretical work treated the 
spin-phonon coupling in terms of a static dimerizationu. 
This approach does not allow the lattice to adjust to spin 
fluctuations and can only be expected to be valid in the 
adiabatic regime. In the case of CuGeOs it has been 
shown that the magnetic energy scales given by the an- 
tiferromagnetic exchange integral J and the phonon fre- 
quencies f2 are of the same order of magnitude^. Thus, 
\-CuGeO3 has to include dynamical 
tB. For technical reasons only very 
little is known about the experimental relevant regime 
with next-nearest, peighbor coupling and intermediate 
phonon couplinglJ'0'E3. In this Letter we present nu- 
merically reliable results in the experimentally relevant 
regime of intermediate phonon frequencies including a 
frustrating next-nearest neighbor interaction. 

We investigate a spin-^ Heisenberg chain coupled to 
dispersionless Einstein phonons with frequency O 



= J 2J(1 + g(b\ + b l ))S l ■ S i+1 + aSi ■ S i+2 ) 



a realistic _mi 
phononsEEll 



+nJ2b\b l (l) 



Here, the b\ and bi are the local phonon creation and an- 
nihilation operators respectively and g is the spin-phonon 
coupling. The ratio of nearest to a frustrating next- 
nearest exchange coupling is given by a. 

There are two physically independent mechanisms in 
this model that can cause a dimerization. First, without 
spin-phonon coupling the frustrating next-nearest neigh- 
bor coupling parameterized by a drives a zero temper- 
ature quantum phase transition from a spin liquid to 
a dimerized phase at a critical value a c . A value of 
a c = 0.241 was determined by numerical studiesEj. In 
the anti-adiabatic limit the phonons can be integrated 
out and the model can be mapped onto a system with 
a frustrating next-nearest neighbor interaction^. For 
a = a phase transition from a gapless spin fluid to 
a gapped dimerized phase occurs at a non-zero value of 
the spin-phonon coupling. These results were confirmed 
to hold in the non-jadiabatic and adiabatic regimes by 
numerical studiesE3'E3o. This in particular is a feature 
of the dynamical model, since the static model exhibits 
a dimerization for all non-vanishing couplingsO. 

Our nu meri cal approach is based on a linked cluster 
expansionlljiHJ which allows to evaluate physical observ- 
ables that are additive when the system separates into 
disconnected parts. It is an inherent feature of cluster 
expansions that physical observables are evaluated in the 
thermodynamic limit. We calculated expansions for the 
ground state energy, Eq, the dispersion of the elementary 
spin triplet excitation Exsi'i) including the spin gap A 
and the structure factor S'(q). The resulting series were 
analyzed using Pade extrapolation techniques. 

In order to perform a dimer expansion of the Hamilto- 
nian (|l|) we divide the Hamiltonian as 
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H = H + H 1 
where we shift the phonon operators 

h = U - —Si ■ S l+1 (2) 

which yields an exactly solvable Hamiltonian Hq 

H = jJ2[Si- s i+1 - f^(3 • s l+1 ) 2 ] + nJ2~ b l~ b i 

i i 
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Here the index i runs over the strong dimer bonds while 
the index j runs over the weak bonds between the dimers. 

The eigenstates of the Hamiltonian Hq are described as 
a product of spin and phonon states. The spin states are 
products of local singlet and triplet dimer states. As a 
consequence of the transformation (|J) the phonon states 
on the strong dimer bonds are coherent states given by 

|0 s/t ) = exp(-r4/2) ■ exp(tj a/t b^)\Q) (4) 

Here r] s — 3<?/4f2 corresponds to a spin singlet state and 
r]t = —g/itt corresponds to a spin triplet state. The weak 
bonds between the dimers are occupied by rij quantum 
phonons, where rij = bjbj. 

The perturbing Hamiltonian Hi is then given by 

H x = AJ^ [(1 + 9{b] + b^Sj ■ S j+1 + aS 3 ■ S j+2 ] (5) 

j 

where the index j runs over the weak bonds between the 
dimers. The expansion parameter A describes the ratio 
of the spin-spin coupling of the weak and strong bonds. 
The expansion is systematic in A while the spin-phonon 
coupling g and the next-nearest neighbor interaction a 
are fixed for each evaluation of a series. 

The total Hilbert-space of the Hamiltonian (|l|) is the 
tensorial product of the space of the spin configurations 
and the phononic space. While there is a finite set of spin 
configurations for every finite cluster, the Hilbert space 
associated with the phonons is infinite even for a finite 
cluster. 

As a starting point for series expansions we use an ini- 
tial wavefunction that describes a dimerized spin phase 
of spin singlets with a macroscopic occupation of phonon 
states on the dimer bonds and no quantum phonons on 
the bonds between the dimers. Performing a series ex- 
pansion, only a limited number of quantum phonons is 
created on the weak bonds between the dimers. The 
macroscopic occupation numbers of the phonon states 
on the strong dimer bonds are tied to the fluctuations of 
the underlying spin states. Thus, the overall phononic 
Hilbert space is finite for a series expansion tn . fi nite or- 
der. In contrast to previous theoretical workllallJ we are 
not forced to truncate the phononic Hilbert space. 



Series expansions were performed up to 10 order with 
ten contributing clusters. The largest cluster contains 
eleven dimers or an equivalent of 22 spins. The Hilbert 
space for this cluster has 9,156,836 contributing states. 
Calculations were performed on local workstations. 

As a closing remark to these technical aspects it is 
stressed again that our calculations do not require any 
finite size scalings. 

The antiferromagnetic exchange integral J and the 
next-nearest neighbor exchange integral J' are deter- 
mined from fits to magnetic susceptibility data. The pa- 
rameter get J = 160 K and a = 0.36 was estimated for 
CuGeOaQ'El. The frequencies of the Peierls-active phonon 
modes were shown to be of the order of the magnetic 
energy scalesa, f2i = 1 J and fl 2 =j— ]2 J. Recently, 
the corresponding coupling constantst3 were proposed 
to be gi — —0.096 and 52 = 0.16. We will investigate 
our model with these parameters. To make contact to 
previous numerical studiesliil a third parameter set with 
0,3 = 0.3 J and 33 = 0.11 is taken into account. 

For the spin-i Heisenberg chain with static dimeriza- 
tion 5 several studies of the gap dependency on the static 
dimerization were performed. The corresponding Hamil- 
tonian is given by 

H> = J o K 1 + (-i)^))^ ■ S i+ i + ao& ■ S i+2 ] (6) 

i 

The parameters of this model are related to those of 
Hamiltonian (]l|) by 

J = J (l + 5) 

A = (1 -$)/(! + *) 

a = a /(l + <5). (7) 

At the critical point the triplet gap A is known to show 
a A(5) tx <5 2 / 3 dependencya. For supercritical frustration 
(a > a c ) a dependency Atf) — A(0) tx <5 2 / 3 was recently 
proposed by Uhrig et. a\E3. 

In order to recover the original Hamiltonian ([!]) the 
evaluated series have to be analyzed at A = 1. We 
have performed two types of Pade approximations. First, 
motivated by the proposed supercritical dependency, we 
evaluate biased Pades fixing the exponent to be 2/3. Sec- 
ond, assuming a power law dependency without restrict, 
ing the exponent, we evaluate Dlog Pade approximantsE3 
and reintegrated the obtained series to evaluate the gap. 

For small supercritical values of a and small coupling 
constants g both Pade approximants show well matching 
results. Varying the degree of numerator and denomina- 
tor polynoms does not substantially change the resulting 
values justifying both approaches. For values of a and 
g far away from critical points the biased Pade approxi- 
mants become less reliable while the Dlog Pade approx- 
imants still give stable results as can be seen in table 
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The inset in Fig. [lj shows the gap in units of Jo versus 
the frustration ctQ. The filled symbols are obtained by 
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FIG. 1. Gap versus spin-phonon coupling. The dashed line 
is a polynomial fit to the data. The inset shows the gap 
versus the frustration in the static model, filled symbols are 
obtained by evaluation of the biased &n5]-Pade approximant. 
The open symbols show DMRG dataN. 
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FIG. 2. Dispersion of the elementary triplet excitation 
evaluated by a dimer expansion up to 10 th order for differ- 
ent phonon frequencies and spin-phonon coupling constants. 
For each momentum the series were analyzed using Dlog Pade 
approximants. 



evaluation of the biased [5, 5]-Eade approximant. The 
open symbols show DMRG dataEj. Our series results are 
in full agreement with those obtained by the perturbative 
approach using flow equationsEj. Good agreement is also 
found with previous DMRG results. 

Fig. [j] shows the gap in units of J for a frustration of 
a = 0.36 for various spin-phonon couplings and phonon 
frequencies. The error bars represent the deviation of 
the accepted Dlog Pade approximants. The dashed lines 
are polynomial fits to the data. For the parameter set 
a = 0.36, f^3 = 0.3 J, and g% = 0.11 we find a triplet gap 
A = (0.172 ± 0.014) J, which is in good agreement with 
previous exact, diagonalization extrapolations for the tt 
phonon modeO. For the parameters relevant to CuGeOa 
our results are below the-experimentally obtained value 
A ps 0.152 J (2.1 meV)EM. 

In the case of a nearly adiabatic phonon frequency 
fl = 0.3 J the gap A monotonically increases with higher 
spin-phonon couplings. For the non-adiabatic phonon 
frequency Q = 1 J the gap shows a gentle downward 
slope. This indicates a suppression of the frustration for 
small phonon frequencies due to a stronger renormaliza- 
tion of the nearest neighbor spin interaction than the 
next-nearest neighbor spin interactions. For higher cou- 
plings this is compensated by an overcritical coupling. 

The dispersion of the elementary triplet excitation is 
shown in Fig. |^. The series expansion method is based on 
Rayleigh-Schrodinger perturbation theory and is there- 
fore supposed to produce well converging results if the 
contributing states are clearly separated by an energy 
gap. In the region |q| > 0.4 tt the series converges very 
fast. For smaller momenta Dlog Pade approximants are 
required to obtain numerically reliable results. 

It turns out that the two parameter sets relevant for 
CuGeOa give nearly matching dispersion relations. The 



maxima of the dispersion relation at q = ±7r are sup- 
pressed in comparison with the known value w(q = 
±7r) = 7r/2 for the unperturbed Heisenberg chain which 
is mainly due to the next-nearest neighbor magnetic in- 
teraction. 

The Peierls ordering structure of the ground state is ac- 
companied by short range magnetic order. Thus, we eval- 
uate the following real-space magnetic correlation func- 
tion 



S 3 = (N Q ■ N 2j+2 ) 



(8) 



where Nj = Sj — Sj+i is the sublattice magnetization 
and Sj denotes the spin on site j. Technically, this is 
done by formally adding a source term 



H 2 = hJ2N 



2j+2 



(9) 



to the Hamiltonian (jlj) and expanding the ground state 
energy to linear order in h. Differentiating with respect 
to h will give the strong coupling expansion for the cor- 
relation function. The determined structure factors are 
well converged. The results presented in the following are 
obtained by direct evaluation of the series without using 
Pade approximants. 
The Fourier transform 



S(q) = ^2e i ^S j 



(10) 



is shown in Fig. p[ Parameters are chosen in the gapped 
phase where the structure factor shows a Lorentzian 
shape which obviously differs from the logarithmic di- 
vergence of the isotropic Heisenberg model at 5(q = 0). 
In the case of weak spin-phonon coupling as given in 
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FIG. 3. Structure factor 5(q) versus momentum q in a 
sixth order calculation. 



TABLE I. [n/m] Dlog Pade approximants to the series for 
the energy gap A for a frustration a — 0.36 and a phonon 
frequency £1 = 1 J. Values are given in units of J. 
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CuGeC>3 this is mostly due to the next-nearest neigh- 
bor interaction. The evaluated series with next-nearest 
neighbor interaction and the series with additional spin- 
phonon coupling cannot be reasonably distinguished. For 
stronger couplings Fig. |^ shows that the spin-phonon in- 
teraction leads to a further diminishing of the maximum 
at q = and a decrease of the spin correlation length. 

In conclusion, series expansion techniques were applied 
to investigate zero temperature properties of a spin-^ 
Heisenberg chain coupled to Einstein phonons which in 
the non-adiabatic limit serves as a model for CuGeC>3. 
We evaluated the dispersion of the elementary triplet ex- 
citation and demonstrated that the spin-phonon coupling 
leads to a renormalization of the dispersion. The de- 
pendency of the triplet gap on the spin-phonon coupling 
shows a non-monotonous behavior in the non-adiabatic 
regime. For the parameters relevant to CuGe03 no sub- 
stantial renormalization of the triplet gap is found. The 
magnetic ordering of the ground state was determined by 
calculation of the spin structure factor. It was found that 
the next-nearest neighbor interaction dominates the spin 
correlations in the case of CuGeC>3. 

We acknowledge helpful discussions on this problem 
with G. S. Uhrig, S. Eggert, R. Werner and C. Gros. 
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